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INTRODUCTION 


Probit  analysis  was  developed  to  analyze  discrete  or  dichotomous  data  obtained 
by  many  research  endeavors  including  natural  or  threshold  response  rate  for  biological 
systems.  We  have  been  using  this  method  as  a  statistical  tool  to  determine  the  probability 
of  dose-response  curves  for  minimum  visible  lesions  (MVLs)  produced  within  the  eye  for 
single  laser  pulses.  In  most  cases,  the  dose  or  laser  pulse  energy  values  required  to 
produce  a  visible  lesion  within  the  eye  are  reported  as  the  ED50  dose  or  that  dose  which 
has  a  probability  of  50%  of  creating  a  visible  lesion.  However,  the  complete  probability 
curve  is  calculated  during  the  iterative  process,  and  the  printout  generally  gives  points 
between  1%  and  99%. 

The  beginning  of  probit  analysis  dates  back  over  60  years  when  the  Karber^ 
method  was  first  published  to  estimate  the  mean  and  standard  deviation  of  the  distribution 
of  data  where  the  distribution  was  unknown.  The  Karber  method  can  still  be  used  as  a 
first  approximation  for  the  nonparametric  estimates  of  the  mean  and  standard  deviation 
for  the  exact  probit  solution  which  involves  a  series  of  successive  approximations.  The 
exact  probit  method  of  solution  has  replaced  most  other  methods  because  it  conveys  the 
maximum  amount  of  information  from  the  data.  In  addition,  the  exact  method  provides 
the  associated  confidence  intervals,  which  allows  inferences  to  be  made  about  the  quality 
of  the  experiment  and  tightness  of  the  data  so  that  obvious  deficiencies  may  be  corrected. 

Graphical  methods^’^  for  the  probit  analysis  have  been  around  for  over  50  years 
and  were  used  quite  extensively  because  of  their  simplicity  and  ease  of  use.  The  data  is 
plotted  using  this  method,  a  straight  line  is  drawn  through  the  data  points,  and  the  ED50 
level  is  read  from  the  graph.  The  slope  of  the  straight  line  can  also  be  determined  from 
the  graph.  Later,  a  simplified  graphical  method  was  developed  which  included  graph 
paper  especially  designed  to  permit  rapid  analysis.  However,  all  of  these  methods  were 
developed  before  the  modem  day  computer  and  therefore  cannot  compete  with  the  speed 
or  accuracy  of  computer-derived  exact  solutions  of  the  probit  iterative  process.  The 
bootstrap  method'^  was  not  considered  in  this  study  because  we  were  interested  in  only 
two-level  trials  (i.e.,  yes/no  data),  and  this  type  of  data  was  not  compared  in  the  study  by 
Foster  and  Bischof.'* 


METHODS 

Most  of  the  original  theoretical  development  of  the  probit  methods  can  be 
attributed  to  Finney^,  who  first  published  a  book  on  probit  analyses  in  1947.  Since  then, 
there  have  been  a  number  of  revisions,  editions,  and  reprinting  of  his  original  book,  and 
in  one  10-year  period  between  1978  - 1988,  there  were  over  2300  citations^  of  his  books 
in  the  literature.  Thus,  his  procedure  is  the  most  widely  used  to  analyze  yes/no  data  and 
other  discrete  event  data.  Finney  developed  the  mathematical  models  to  utilize  probit 
analyses,  but  it  still  takes  a  computer  program  to  process  the  data.  Many  programs  for 
many  different  computers  have  been  written  to  perform  the  probit  procedures  which  were 
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developed  by  Finney,  and  most  will  calculate  the  probability  curve  equally  well. 

However,  to  standardize  the  reporting  of  data  in  the  literature,  the  S  AS  Institute,  Inc.^, 
Cary,  NC,  developed  computer  programs  for  many  different  computer  systems  and 
operating  systems  which  calculate  maximum-likelihood  estimates  of  regression 
parameters  and  natural  (threshold)  response  rate  for  the  discrete  event  data.  Thus,  these 
programs  calculate  the  exact  same  results  for  all  systems.  One  of  the  authors  of  this 
report.  Captain  Lonnie  Manning,  has  written  a  program  in  C++  that  operates  on  a  PC 
computer  to  utilize  the  methods  developed  by  Finney  (see  Appendix  B).  This  program 
produces  all  of  the  same  output  parameters  as  SAS  Probit.  The  SAS  Institute  has  named 
its  process  The  Probit  Procedure,  and  we  call  ours  EZ-probit. 

There  were  several  reasons  for  writing  the  EZ-probit  program.  In  the  literature, 
we  find  many  references  to  Finney  and  the  method  he  employs  for  the  Probit  calculations. 
However,  SAS's  Probit  and  Finney's  Probit  are  different.  Although  Finney  does  briefly 
mention  the  method  employed  by  SAS  (SAS's  Probit  is  actually  a  Normit  since  the  Probit 
is  not  centered  about  Y-5).  The  EZ-probit  program  incorporates  Finney's  expressions 
exactly.  Comparison  of  the  SAS  output  and  that  from  Finney's  Probit  shows  they  differ 
for  a  small  number  of  data  points  wiA  agreement  increasing  as  the  number  of  data  points 
increase.  This  is  most  noticeable  in  the  fiducial  limit  calculations.  In  most  instances,  the 
difference  is  inconsequential  compared  to  the  experimental  errors  inherent  in  our  type  of 
experiments(i.e.+/-  7%  for  energy). 

The  Probit  Procedure  computes  maximum-likelihood  estimates  of  the  slope  and 
intercept  of  the  probit  equation  using  a  modified  Newton-Raphson  algorithm.  The  data 
set  used  by  SAS  Probit  must  include  either  a  response  variable  giving  the  level  of 
response  for  each  observation  or  a  pair  of  variables  giving  the  number  of  subjects  tested 
and  the  number  of  subjects  responding  for  each  dose  of  the  independent  variable  values. 
Two  goodness-of-fit  cld-square  values  are  computed  if  requested.  Inverse  confidence 
limits  for  one  of  the  independent  variables  can  be  requested,  and  the  confidence  limits  are 
computed  using  a  critical  value  of  1.96,  which  corresponds  to  an  approximate  95% 
confidence  interval.  If  the  Pearson  goodness-of-fit  chi-square  test  is  requested  and  the  p- 
value  for  the  test  is  too  small,  variances  and  covariances  are  adjusted  by  a  heterogeneity 
factor,  and  a  critical  value  from  the  t-distribution  is  used  to  compute  the  fiducial  limits 
(FLs).  The  p-value  used  for  the  chi-square  test  can  be  set  to  different  levels  with  a  default 
p-value  of  0.10.  Also  calculated  and  outputted  is  the  slope  of  the  probit  line  between  the 
ED84  and  the  ED50  values. 

The  EZ-probit  program  is  faster  and  much  more  user-friendly  than  the  existing 
SAS  program.  All  relevant  parameters  can  be  set  by  the  click  of  a  mouse.  To  use  SAS's 
Probit,  the  whole  SAS  program  must  be  loaded.  EZ-probit  is  a  Windows  application 
written  in  C++  (Borland  4.0  Compiler)  and  comes  in  both  a  16  bit  or  32  bit  application. 

In  addition,  EZ-probit  allows  flexibility  in  setting  the  confidence  to  something  other  than 
95%  which  is  currently  very  difficult  to  accomplish  in  SAS. 
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A  good  treatise  on  the  graphical  methods  was  developed  by  Frisch.^  Step-by-step 
procedures  are  presented  in  that  report,  together  with  comparisons  between  graphical, 
Karber,  and  exact  probit  solution.  In  this  paper,  we  present  his  results  and  compare  them 
to  results  run  on  S  AS  Probit  and  EZ-Probit,  along  with  results  from  the  MVL  and  other 
data  sets  created  to  influence  their  ou^uts. 


DATASETS 

The  first  data  set  presented  comes  from  Frisch^  which  is  for  the  70-millisecond 
(ms)  pulse  duration  and  170  exposures.  In  his  report,  Frisch  thoroughly  analyzes  the 
three  methods  and  clearly  explains  each  step  in  each  method.  Real  data  are  used  in  each 
case.  His  report  should  be  consulted  for  learning  each  procedure  and  what  auxiliary  data 
is  necessary.  Table  1  lists  the  results  from  Frisch’s  calculations  for  the  graphical  method, 
Karber  approximation,  his  exact  probit  solution,  and  die  results  from  SAS  Probit  and  our 
EZ-Probit  (both  of  which  calculate  the  exact  probit  solution  using  an  iterative  process). 

Table  1 .  Comparison  of  Four  Methods  for  Data  from  Frisch.^ 


Parameter 

Graphical 

Karber 

Exact  Probit  SAS  Probit 

EZ-Probit 

EDso 

12.5 

13.15 

12.5 

12.4 

12.4 

Upper  FL 

14.4 

14.55 

14.0 

13.8 

13.8 

Lower  FL 

11.2 

11.65 

11.0 

10.8 

10.8 

Slope  of  Probit 

5.9 

— 

5.68 

5.67 

5.67 

Ratio(ED84  /  ED50)  1 .48 

.... 

1.50 

1.50 

1.50 

Analysis  of  Table  1 

The  listed  graphical,  Karber,  and  exact  probit  parameter  values  above  were 
calculations  by  Frisch,  and  the  SAS  Probit  was  run  on  the  PC  version  of  The  Probit 
Procedure.  The  exact  probit  used  by  Frisch  should  be  equivalent  to  the  SAS  probit  since 
they  both  use  the  iterative  process  to  achieve  a  final  solution.  The  differences  between 
the  three  methods  (graphical,  Karber,  and  exact  probit)  are  minimal,  and  it  is  obvious  that 
the  three  methods  give  comparable  results.  SAS  Probit  and  the  EZ-Probit  also  print  out 
the  goodness-of-fit  test,  and  for  the  Pearson  chi-square,  the  value  was  2.23  and  the 
Prob>Chi-Sq  was  0.97.  SAS  Probit  and  EZ-Probit  gave  identical  results  out  to  four 
decimal  places.  Therefore,  this  data  set  did  have  a  good  distribution.  However,  this  is 
only  one  data  set,  and  the  results  of  the  calculations  depend  on  the  data  set.  One 
difference  between  Frisch  and  the  SAS  Probit  is  the  way  slope  is  defined  and  printed  out 
by  SAS  Probit.  Frisch  defines  slope  of  the  probit  as  the  ratio  of  ED84  /  ED50,  while  SAS 
Probit  defines  it  as  the  slope  of  the  straight  line  of  best  fit  to  the  data.  The  two  are 
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inversely  related,  and  the  slope  of  the  probit  may  be  obtained  from  the  ED84  /  ED50  ratio 
simply  by  taking  the  reciprocal  of  the  logariflimio  of  the  ratio  (ED84  /  ED50).  About  2.5  is 
the  crossover  point  of  these  two  numbers,  and  the  theoretical  minimum  for  the  ratio  is  “1” 
(i.e.,  the  slope  of  the  probit  would  be  infinite  in  this  case).  However,  for  made-up  data, 
the  smallest  ratio  is  1.02,  while  the  slope  is  calculated  to  be  90  (see  Table  9). 

In  order  to  make  a  valid  comparison,  real  data  sets  as  well  as  created  data  sets 
must  be  analyzed  using  all  of  the  different  methods,  and  their  outputs  must  be  compared. 
Other  data  sets  compared  were  the  MVL  data  for  90-  and  600-femtosecond  (fs) 
pulsewidths.  These  data  sets  consisted  of  122  exposures  at  90  fs  and  121  exposures  at 
600  fs,  both  at  580-nanometer  (nm)  wavelength  for  the  rhesus  monkey  eye.  A  graphical 
analysis  was  performed  on  each  set  and  compared  with  SAS  Probit  results.  The  results 
for  Aese  two  sets  are  listed  in  Tables  2  and  3. 

Table  2.  MVL  Data  at  90  Femtoseconds. 


Parameter 

Graphical 

Karber 

SAS  Probit 

EZ-Probit 

ED50 

0.46 

0.60 

0.43 

0.43 

Upper  FL 

0.77 

0.71 

0.61 

0.61 

Lower  FL 

0.25 

0.49 

0.27 

0.27 

Slope  of  Probit 

1.04 

— 

1.61 

1.61 

Ratio(ED84  /  ED50) 

9.9 

4.2 

4.2 

Analysis  of  Table  2 

Table  2  again  shows  the  comparison  between  graphical,  Karber,  and  exact  probit 
as  calculated  by  SAS  Probit  and  EZ-Probit.  This  data  set  was  selected  because  of  its 
spread  in  data  and  its  low  slope  calculated  by  SAS.  This  data  set  illustrates  the  necessity 
of  the  iterative  process  involving  a  series  of  successive  approximations  to  obtain  the  final 
solution.  The  Karber  method  gives  only  a  first  approximation,  while  die  graphical 
method  gives  only  an  eyeball  approximation  with  its  linear  regression.  SAS  Probit 
calculated  a  higher  slope  value  than  the  graphical,  1.61  versus  1.04,  respectively.  The 
goodness-of-fit  test  gave  a  value  of  55  and  a  Prob>Chi-Sq  of  0.92.  If  one  assumes  that  a 
good  and  reasonable  data  set  should  have  at  least  a  0.90  for  a  Prob,  then  this  data  set  does 
fit  the  requirement. 

The  data  set  shown  in  Table  3  was  selected  because  of  its  high  value  of  slope  of 
the  probit  line  and  tightness  of  data  points.  The  goodness-of-fit  test  gave  a  value  of  30 
with  a  Prob>Chi-Sq  of  1.0000.  In  this  case,  the  three  methods  gave  very  close  results 
which  is  a  good  indicator  of  the  power  of  the  interactive  process  of  SAS  as  compared  to 
the  data  in  Table  2.  The  other  two  methods  give  good  agreement  in  their  approximation 
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because  of  the  data  set.  Not  only  are  the  EDsos  close  between  the  three  methods,  but  the 
FLs  are  also  very  close  to  each  other.  Thus,  any  of  these  methods  would  suffice  if  all  data 
sets  were  as  regular  as  this  set. 

Table  3.  MVL  Data  for  600  Femtoseconds. 


Parameter 

Graphical 

Karber 

SAS  Probit 

EZ-Probit 

ED50 

0.24 

0.28 

0.26 

0.26 

Upper  FL 

0.32 

0.32 

0.31 

0.31 

Lower  FL 

0.18 

0.24 

0.21 

0.21 

Slope  of  Probit 

3.31 

— - 

4.79 

4.79 

Ratio(ED84  /  ED50) 

2.40 

1.58 

1.58 

Analysis  of  Table  3 

Some  questions  have  arisen  concerning  whether  the  distribution  of  data  in  the  set 
is  normal  or  unknown.  Since  the  Karber  method  reportedly  does  not  depend  on  a  known 
distribution,  we  attempted  to  show  the  differences  between  the  methods  by  using  several 
known  distributions  of  data.  Thus,  chi-squared,  linear,  step,  and  random  distributions 
were  used  in  the  following  data  sets. 

The  data  presented  in  Tables  4  and  5  were  obtained  from  the  equation  Y  =  5.6  + 
2x,  which  plots  as  a  straight  line  on  linear  graph  paper.  The  dosage  or  trial  level  point 
distribution  was  obtained  from  a  table  of  chi-square  values,  and  fliere  is  no  spread  in  the 
data.  This  normal  chi-square  distribution  of  90  doses  or  levels  with  linear  increasing  yes 
responses  is  presented  in  Table  4,  and  the  same  distribution  with  varying  number  of  data 
points  is  presented  in  Table  5. 

Table  4.  Comparison  of  Three  Methods  for  an  Exact  Straight  Line  Distribution. 

(Ninety  doses  with  chi-square  distribution  of  levels.) 


Parameter  Graphical  Karber  SAS  Probit  EZ-Probit 


ED50 

0.50 

0.50 

0.50 

0.50 

Upper  FL 

0.70 

0.72 

0.72 

0.72 

Lower  FL 

0.36 

0.40 

0.35 

0.35 

Slope  of  Probit 

2.0 

— 

2.0 

2.0 

Ratio(ED84  /  ED50) 

3.2 

— - 

3.2 

3.2 

5 


Analysis  of  Table  4 


This  data  set  was  chosen  to  see  how  well  SAS  Probit  performed  as  compared  to 
the  graphical  procedure  since  this  is  a  perfect  data  set  with  chi-square  distribution  and  no 
spread.  In  this  case,  the  graphical  procedure  should  give  an  exact  solution  and  the  KSrber 
method  could  also  be  compared.  The  goodness-of-fit  test  gave  a  value  of  0.0021  with  a 
Prob  of  1.0000.  As  can  been  seen  from  the  results  in  the  table,  there  is  very  good 
agreement  between  the  four  methods.  Thus,  for  this  type  of  data  all  three  methods 
perform  equally  well. 

Table  5.  Same  Distribution  as  Table  4  for  Different  Number  of  Data  Points. 

(All  calculations  with  SAS  Probit.  EZ-Probit  values  shown  in  ( ).) 


Parameter  ( #  nts) 

180 

20 

45 

22 

18 

14 

ED50 

0.5 

0.5 

0.49 

0.53 

0.50 

0.50 

Upper  FL 

0.64 

0.72 

0.74 

1.14 

1.45 

(1.32) 

9.54 

(1.67) 

Lower  FL 

0.39 

0.35 

0.32 

0.27 

0.17 

(.19) 

0.02 

(.15) 

Slope  of  Probit 

2.0 

2.0 

2.8 

2.4 

2.8 

5.3 

Ratio(EDg4  /  ED50) 

3.3 

3.3 

2.3 

2.9 

2.4 

1.44 

HD99 

7.31 

7.31 

3.35 

4.92 

3.43 

1.38 

EDoi 

0.03 

0.03 

0.07 

0.06 

0.07 

0.18 

G-o-F  Value 

0.0042 

0.0021 

1.25 

3.39 

2.67 

1.6 

Prob>Chi-Sq 

1.00 

1.00 

0.99 

0.85 

0.91 

0.98 

Analysis  of  Table  5 

The  data  sets  presented  in  Table  5  were  all  run  with  SAS  Probit  and  EZ-Probit  to 
show  the  effects  of  sample  size  on  the  various  calculated  parameters.  EZ-Probit  gave 
identical  results  to  SAS  Probit  with  the  exceptions  of  the  two  values  for  the  fiducial  limits 
shown  in  parentheses  below  each  value.  Since  the  EZ-probit  program  incorporates 
Finney's  expressions  exactly,  there  is  slight  differences  between  the  two  methods  for 
calculating  the  fiducial  limits  as  shown  by  these  small  differences.  In  most  instances 
however,  the  differences  are  inconsequential  compared  to  the  experimental  errors 
inherent  in  our  type  of  experiments(i.e.+/-  7%  for  energy).  The  variations  in  the  ED50 
values  were  due  to  the  inability  to  maintain  the  exact  percentages  of  yes  responses  in  each 
level  because  of  so  few  points  at  each  level.  With  the  fewer  points  (below  45),  the 
percentages  at  each  level  were  no  longer  constant  and  there  was  spread  in  the  data  set  so 
that  not  all  points  fell  on  the  straight  line.  From  this  table,  it  can  be  seen  that  the 
midpoint,  ED50 ,  remains  fairly  constant  while  the  FLs  spread  apart  as  the  number  of 
points  decreases.  In  fact,  below  14  data  points,  neither  program  could  obtain  FLs  with  a 
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95%  confidence  level.  However,  it  is  possible  to  run  EZ-Probit  at  any  confidence  level 
specified.  It  is  clear  that  the  slope  of  the  probit  changed  because  of  the  fewer  number  of 
“T’s  below  0.5  and  the  ones  used  were  closer  to  the  ED50  point  (i.e.,  smaller  spread  in  the 
data).  The  span  between  the  EDoi  and  ED99  points  decreased  as  the  number  of  points 
decreased  because  of  the  smaller  spread  in  data.  From  these  data  sets,  it  appears  that 
“reasonable”  FLs  are  not  obtained  until  the  data  set  is  larger  than  27.  Herein  we  define 
reasonable  to  be  when  the  upper  FL  is  no  larger  than  50%  of  ED50,  and  the  lower  FL  is  no 
smaller  than  50%  of  ED50. 


The  next  data  set  (Table  6)  to  be  compared  by  the  three  methods  was  a  set  of  ten 
dosage  levels,  equally  spaced,  with  a  linear  increase  from  0  to  100%  yes  responses  for  the 
levels.  The  levels  increased  from  0.1  to  1.0  in  0.1  increments,  and  the  response  increased 
proportionally. 

Table  6.  Comparison  for  Equal  Doses,  Equally  Spaced  and  Linear  Response  Levels  by 
Four  Probit  Methods  for  95  Trials.  SAS  Probit  and  EZ-Probit  equal  to  four 
decimals. 


Parameter 

Graphical 

Karber 

SAS  Probit 

EZ-Probit 

ED50 

0.46 

0.50 

0.50 

0.50 

Upper  FL 

0.63 

0.56 

0.60 

0.60 

Lower  FL 

0.40 

0.44 

0.41 

0.41 

Slope  of  Probit 

4.0 

4.0 

4.0 

Ratio(EDg4  /  ED50) 

1.81 

.... 

1.91 

1.91 

ED99 

1.91 

1.91 

EDoi 

0.13 

0.13 

G-o-F  Value 

2.13 

2.13 

Prob>Chi-Sq 

0.98 

0.98 

Analysis  of  Table  6 


By  comparing  these  four  results,  all  methods  work  equally  weU  on  this  data  set. 
However,  a  straight  line  was  drawn  through  a  severely  curved  set  of  data  points,  and  it 
was  not  obvious  by  looking  at  the  graph  that  it  would  give  correct  results.  In  this  case, 
eyeballing  gave  results  as  good  as  any  of  the  other  three  methods.  Thus,  equally  spaced 
dosages  or  trial  levels  appear  to  be  as  good  as  or  better  than  chi-square  distribution  of 
levels.  The  goodness-of-fit  test  gave  a  value  of  2.13,  and  the  Prob  was  0.98  for  both 
methods  run  on  the  computer. 

The  next  set  of  calculations  was  performed  on  the  same  distribution  but  with  a 
decreasing  number  of  samples.  The  same  percentages  of  yes  responses  were  maintained 
to  their  limit.  These  calculations  were  all  performed  with  SAS  Probit  and  EZ-Probit 
(Table  7),  and  the  sample  size  was  reduced  from  95  data  points  to  12.  EZ-Probit  values. 
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if  different  from  SAS  Probit  values,  are  shown  in  parentheses.  Note  that  where  SAS 
Probit  was  not  able  to  calculate  the  FLs  at  the  95%  confidence  level  for  the  12  data 
points,  EZ-Probit  was  able  to  calculate  the  limits  at  the  85%  confidence  level. 


Table  7.  A  Comparison  of  SAS  Probit  Calculations  for  Data  Sets  with  Decreasing 
Number  of  Data  Points.  (Same  data  as  in  Table  6.) 


Parameter  (data  pts)  95  45  22  .14  13  12 


ED50 

0.50 

0.49 

Upper  FL  0.60 

0.62 

0.65 

Lower  FL 

0.41 

0.37 

Slope  of  Probit 

4.0 

5.1 

Ratio(ED84  /  ED50) 

1.81 

1.57 

ED99 

1.91 

1.40 

EDoi 

0.13 

0.17 

G-o-F  test  Value 

2.13 

1.46 

Prob>Chi-Sq 

0.98 

0.98 

0.49 

0.49 

0.47 

0.49 

0.76 

0.76 

none 

(.64) 

(.72) 

(.72) 

(.63@85%) 

0.36 

0.20 

0.002 

none 

.36) 

(.24) 

(.10) 

(.34@85%) 

10.2 

8.7 

8.0 

7.8 

1.27 

1.28 

1.36 

1.35 

0.84 

0.90 

0.93 

0.97 

0.29 

0.26 

0.24 

0.24 

0.12 

0.73 

1.17 

0.48 

0.99 

0.98 

0.95 

0.99 

Analysis  of  Table  7 

From  this  data  set,  it  appears  that  the  minimum  number  of  data  points  to  have 
reasonable  FLs  for  this  distribution  would  be  about  22  data  points  because  even  by 
doubling  or  quadrupling  this  number,  the  FLs  change  minimally.  With  only  12  data 
points,  SAS  could  not  calculate  the  FLs,  and  the  limits  were  extremely  wide  with  only  13 
points.  With  14  points,  however,  the  FLs  were  stiU  not  reasonable,  but  were  very  close  to 
the  limits.  Again,  there  are  slight  differences  between  the  fiducial  Hmits  at  the  lower 
number  of  data  points  which  can  be  attributed  to  the  different  methods  in  calculating  the 
fiducial  limits. 

Next,  data  points  were  added  to  the  14-pomt  set  in  an  attempt  to  reduce  the  spread 
in  the  FLs  to  equal  those  for  the  95-point  data  set.  Table  8  gives  the  results  of  adding  8 
data  points  in  various  combinations  to  the  14-point  data  set  in  Table  7  above.  The  desire 
was  to  minimize  the  total  number  of  exposures  or  samples  and  yet  obtain  the  large 
sampling  results  as  those  obtained  for  95  data  points.  The  starting  data  was  as  follows: 

dose  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

shots  1  2  3  2  3  2  1 

#  yes  0  0  1  1  2  2  1 
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Table  8.  Additional  Data  Points  Added  to  Optimize  the  FLs  and  Minimize  the 
Number  of  Trials.  (Add  eight  data  points:  four  data  points  added  at  ED25  and 
four  data  points  added  at  ED75  or  at  EDos  and  at  ED92.) 


Data  pts  added 

EDjn 

Upper  FL 

Lower  FL 

Slope 

G-of-F 

Chi-Sq 

Original  95  pts 

0.50 

0.60 

0.41 

4.0 

2.13 

0.98 

Original  22  pts 

0.49 

0.65 

0.36 

10.2 

0.12 

0.99 

Original  14  pts 

0.49 

0.76 

0.20 

8.7 

0.73 

0.98 

4-0.5@‘T” 

4-0.5(®“0” 

0.49 

0.69 

0.27 

8.7 

0.75 

0.98 

4 -0.4  @75%  -’T” 

4 -0.6  @25%  -”1” 

0.48 

none  none  4.0 

(0.71@85%)  (0.31@85%) 

4.16 

0.52 

4 -0.4  @50%  -’T” 

4 -0.6  @50%  -”1” 

0.48 

0.79 

0.25 

5.7 

1.83 

0.87 

4 -0.4  @25%  -”1” 

4 -0.6  @75%  -’T” 

0.49 

0.62 

0.36 

8.3 

0.57 

0.99 

4-0.4@0%  -”1” 

4 -0.6  @100%-”!” 

0.49 

0.58 

0.40 

12.7 

0.09 

1.00 

4 -0.3  @25%  -’T” 

4 -0.7  @75%  -”1” 

0.47 

0.67 

0.32 

5.5 

0.24 

1.00 

4-0.3@0%  -”1” 
4- 0.7  @100%-”!” 

0.48 

0.59 

0.36 

10.5 

1.44 

0.92 

Analysis  of  Table  8 

With  95  data  points,  the  FLs  were  0.60  and  0.41  as  compared  to  the  0.76  and  0.20 
with  only  14  data  points.  When  all  eight  data  points  were  added  at  the  ED50  level,  four 
‘T”s  and  four  “0”s,  the  FLs  changed  very  little.  Next,  four  data  points  were  added  at  the 
ED25  level  and  four  data  points  added  at  the  ED75  level  in  various  combinations  as  shown 
in  the  table.  With  three  ‘T”s  and  a  zero  added  at  0.4  level,  and  three  “0”s  and  a  “1”  added 
at  0.6  level,  the  FLs  could  not  be  calculated  with  S  AS  because  the  probit  curve  had 
predicted  the  reverse  probabilities  (i.e.,  three  “0”s  and  a  “1”  at  0.4,  and  the  reverse  at  0.6). 
EZ-Probit  was  able  to  calculate  the  FLs  at  the  85%  confidence  level,  and  the  limits  are 
shown  in  parentheses.  It  should  be  pointed  out  that  the  Prob>Chi-Sq  was  only  0.52, 
which  showed  very  poor  correlation.  Next,  two  ‘T”s  and  two  “0”s  were  added  at  each 
point,  which  was  still  above  the  predicted  values  and,  therefore,  the  FLs  were  not  helped. 
Also,  the  Prob>Chi-Sq  was  still  less  than  0.90,  which  indicated  that  more  data  would  be 
needed  before  being  acceptable.  Whenever  the  predicted  results  were  added  (i.e.,  three 
“0”s  and  one  “1”  added  at  0.4,  and  three  ‘T”s  and  one  “0”  added  at  0.6,  the  FLs  were 
drawn  closer  together  and  the  Prob  became  0.99).  When  four  “0”s  were  added  to  0.4  and 
four  ‘T”s  were  added  to  0.6  levels,  the  FLs  were  actually  closer  fliat  those  obtained  with 
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95  data  points.  Other  data  points  in  two  different  combinations  were  added  to  the  EDog 
(level=0.3)  and  an  ED92  (level=0.7)  predicted  levels  as  shown  in  die  table.  In  both  cases, 
the  FLs  were  brought  closer  together  but  not  as  much  as  before. 

CAUTION:  It  must  be  pointed  out  that  the  additional  data  added  to  the  above  table  was 
always  added  with  equal  number  of  data  points  above  and  below  the  ED50  values  at 
dosage  level  exactly  spaced  above  and  below  the  ED50  value.  This  symmetry  of  adding 
data  is  necessary  to  prevent  having  a  biased  estimate  of  the  ED50.  Thus,  data  may  not  be 
arbitrarily  added  to  an  existing  data  set  without  carefully  designing  the  additional 
experimental  data  sets. 


Another  data  set  which  is  encountered  sometimes  with  real  data  is  the  step 
function  response  in  which  the  responses  (0  or  1)  are  all  “0”  below  a  certain  level  or  dose 
and  all  “1”  above  the  level  or  dose.  This  data  set  cannot  be  analyzed  graphically  or  by  the 
Karber  method  because  the  slope  of  the  probit  curve  would  be  a  vertical  straight  line,  and 
no  points  would  fall  on  the  graph.  Also  included  was  an  increasing  number  of  points 
overlapping  from  1  and  2  points  to  a  constant  or  flat  response  of  50%  between  two  levels. 
In  Table  9,  the  response  of  p  =  .5  is  constant  for  levels  of  0.4, 0.5,  and  0.6  in  the  4th 
column.  Columns  5  and  6  both  have  p  =  .5  for  levels  .4,  .5,  and  .6  and,  in  addition,  levels 
.3  and  .7  have  either  3  out  of  10  shots  or  5  out  of  10  shots  for  p  =  .3  or  p  =  .5.  The  data  in 
column  5  with  p  =.3  for  level  0.3  is  30%  (i.e.,  three  ‘T”s  out  of  10  trials  and  three  “0”s  in 
level  .7,  while  in  the  5th  column  at  p  =  .5,  the  response  is  50%  for  both  the  0.3  and  0.7 
levels).  In  other  words,  the  flat  response  occurs  between  0.3  and  0.7  levels  or  for  5  levels 
in  column  6.  All  computer  runs  with  SAS  contained  100  data  points. 

Table  9.  SAS  Probit  Calculations  for  Step  Functions  and  Constant  Response  Levels. 
Step  distribution  at  0.5  and  p  =  .5  at  .3,  .4,  .5,  .6,  and  .7.  One  hundred  data 
points  used.  SAS  Probit  and  EZ-Probit  gave  identical  results  as  shown. 


Parameter 

Step  function 
no  overlap 

1  pt  2  pts 
overlap  overlap 

p=.5 

t.4..5..6) 

p=.3 

p=.5 

ED50 

0.5 

0.51 

0.49 

0.48 

0.46 

0.44 

Upper  FL 

none 

0.55 

0.54 

0.53 

0.53 

0.59 

Lower  FL 

none 

0.46 

0.45 

0.42 

0.38 

0.28 

Slope  of  Probit 

89.8 

20.6 

15.4 

8.2 

5.2 

4.2 

Ratio(ED84  /  ED50) 

1.02 

1.12 

1.18 

1.33 

1.56 

1.7 

ED99 

0.53 

0.66 

0.70 

0.92 

1.28 

1.60 

EDoi 

0.47 

0.39 

0.35 

0.25 

0.16 

0.12 

G-o-F  Value 

0.00 

0.47 

0.21 

10.1 

8.46 

16.2 

Prob>Chi-Sq 

1.00 

1.00 

1.00 

0.26 

0.39 

0.04 
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Analysis  of  Table  9 


With  no  overlapping  of  data  points,  neither  SAS  Probit  nor  EZ-Probit  could 
calculate  FLs,  and  the  calculated  slopes  from  both  were  very  large.  Also  note  that  the 
ED99  and  EDoi  are  almost  the  same  as  the  ED50.  As  the  number  of  overlapping  data 
points  increases  and  the  slope  of  the  probit  decreases  very  rapidly,  there  is  not  much 
change  in  the  fiducial  limits.  For  the  last  two  columns,  adding  “P’s  to  the  0.3  level 
decreased  the  ED50  value  and  increased  the  spread  between  the  FLs.  This  occurs  because 
the  number  of  ‘T”s  and  “0”s  added  did  not  match  predicted  values  for  the  levels  .3  and  .7. 
Also  note  that  in  the  last  two  columns  the  value  was  greater  than  10  and  the  Prob  was 
between  0.39  and  0.04.  These  are  very  different  values  from  those  calculated  for  a 
normal  distribution  and  would  indicate  that  more  data  points  need  to  be  taken  even 
though  the  FLs  may  appear  to  be  reasonable. 


CONCLUSIONS 

The  main  conclusion  to  be  drawn  from  this  study  is  that  the  exact  probit  procedure 
is  the  one  that  should  be  used  on  any  set  of  data  and  that  the  graphical  solution  is  not 
adequate  with  todays  computer  technology.  As  quoted  by  Frisch^  “Graphical  solutions 
very  often  are  adequate  but  complications  may  make  this  type  of  solution  unsatisfactory  — 
and--  the  exact  probit  method  of  solution  is  the  most  advantageous  since  it  conveys  the 
maximum  amount  of  information  from  the  data”.  Thus  with  the  computer  solutions 
available  today,  it  is  unthinkable  to  resort  to  graphical  or  calculator  solutions  such  as  the 
Karber  method  for  only  a  first  approximation. 

Many  other  conclusions  may  be  drawn  from  these  data  sets,  but  the  main 
consideration  is  the  ability  to  maximize  the  reliability  of  the  experimental  procedure  with 
a  minimum  number  of  data  points.  It  was  shown  that  the  ED50  may  be  obtained  with  FLs 
at  their  95%  confidence  interval  with  as  few  as  13  data  points,  even  though  the  span 
between  the  FLs  may  be  wide.  It  was  also  shown  that  with  only  eight  more  data  points 
added,  the  FLs  could  be  as  good  as  or  better  than  if  95  data  points  had  been  taken, 
depending  on  the  data  points  added.  It  may  be  concluded  that  when  data  points  are  added 
to  an  existing  data  set  with  95%  FLs,  whenever  those  data  points  are  at  doses  below  the 
ED50  values  and  the  percentages  of  yes  or  “1”  is  less  than  the  ED  percentage  at  that  dose, 
the  combined  data  set  will  be  helped.  If  the  percentages  of  yes  or  “1”  are  greater  than  the 
ED  value,  the  data  set  will  not  be  helped.  In  other  words,  for  a  fixed  number  of  data 
points  added  below  the  ED50,  the  percentage  of  positive  responses  must  be  less  than  the 
ED  percent  at  that  dose  or  the  lower  FL  will  be  further  from  the  ED50  value  and  may  not 
be  calculable  at  the  95%  confidence  level.  The  reverse  is  hue  for  data  points  added  above 
the  ED50  where  the  percentages  of  yes  or  “1”  must  be  greater  than  the  ED  percentage  at 
that  dose  level.  The  reasoning  for  the  above  is  obvious  due  to  the  maximum  likelihood 
estimators  of  the  probit  procedure,  and  when  the  new  data  does  not  match  the  probit 
probabilities,  the  combined  data  set  becomes  either  better  or  worse. 
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Another  conclusion  with  the  most  significance  is  that  it  is  relatively  easy  to  write 
your  own  computer  program  to  calculate  the  probability  curve  and  the  FLs  for  any  level 
of  significance  desired.  EZ-Probit  has  been  shown  to  be  equivalent  to  the  SAS  Probit 
with  the  exception  of  the  fiducial  limits  for  a  limit  number  of  data  points,  and  also  has  the 
added  flexibility  to  be  able  to  determine  the  FLs  at  any  desired  level  of  confidence.  EZ- 
Probit  is  a  lot  easier  to  use,  faster  than  SAS  Probit,  and  has  more  capability.  Since  EZ- 
Probit  has  been  installed  on  the  P-drive,  it  is  available  for  use  by  anyone,  anywhere, 
anytime. 
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APPENDIX  A: 

PC  COMPUTER  PROGRAM 
FOR 

SAS  PROBIT 
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The  following  program  was  written  for  SAS  to  operate  on  the  PC  computer  and  prints  the 
outputs.  A  brief  explanation  of  what  the  code  does  follows  the  program.  A  complete 
description  of  the  terminology  is  given  in  the  SAS/STAT  Software  Syntax  manual.^ 

1.  filename  testdata  ‘filename.dat’; 

2.  title  'filename.dat'; 

3.  data  indata; 

4.  infile  testdata; 

5.  input  energy  mvl; 

6.  run; 

7. 

8.  proc  sort  data  =  indata; 

9.  by  energy; 

10.  run; 

11. 

1 2.  proc  print  data  =  indata; 

13. 

14.  run; 

15. 

1 6.  proc  sort  data  =  indata; 

17.  by  desending  energy; 

18.  run; 

19. 

20.  proc  probit  data  =  indata  outest  =  mvldata  order  =  data  covout  Hprob  =  0. 1 0 
lackfit  log  10  inversecl; 

21.  class  mvl; 

22.  model  mvl  =  energy  /  lackfit  d  =  normal  corrb  covb  inversecl; 

23.  output  out  =  mdat  prob  =  p  xbeta  =  xB  std  =  sd; 

24.  run; 

25. 

26.  goptions  target=wmprtg;  /*  grey-scale  postscript  */ 

27. 

28.  goptions  rotate  =  landscape; 

29. 

30.  proc  gplot; 

3 1 .  axis  logbase  =10  logstyle  =  power; 

32.  plot  mvl*energy  p*energy  /  overlay; 

33.  run; 


The  first  6  lines  read  the  data  from  a  file  into  a  SAS  data  set  variable  named  “indata.” 
The  file  should  be  a  two-column  list  of  numbers  separated  by  a  space  with  each  line 
terminated  with  a  return  character.  The  first  colunm  should  be  the  dose  or  energy 
delivered  and  will  be  read  into  the  “energy”  array.  The  second  column  should  be  the 
yes/no  response  (1  =  yes,  0  =  no)  to  the  dose.  This  colunm  will  be  read  into  the  “MVL” 
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array.  Lines  8  through  10  sort  the  data  set  by  energy  in  ascending  order  for  printout  in 
lines  12  through  14.  Lines  16  through  18  sort  the  data  again  by  energy  in  descending 
order  for  the  probit  procedure.  This  last  sort  is  to  ensure  that  the  first  data  point  is  a  yes 
response.  The  SAS  probit  procedure  assumes  that  the  first  data  point’s  response  is  the 
value  that  will  be  considered  a  “yes”  response  for  the  current  run.  The  call  to  the  probit 
procedure  begins  on  line  20  and  ends  on  line  24.  The  following  is  a  description  of  the 
flags  used  in  the  call  to  the  probit  procedure  and  their  purpose. 


data 

outest 

covout 

Hprob 

lackfit 

loglO 

inversecl 


Supplies  the  name  of  die  SAS  data  set  to  be  used  by  the  probit  procedure. 
Supplies  the  name  of  the  SAS  data  set  where  the  output  will  be  written. 
Writes  the  covariance  estimates  to  the  outest  data  set. 

Specifies  the  probability  level  to  indicate  the  goodness  of  fit. 

Performs  a  Pearson  chi-square  test  of  the  fit  and  log-likelihood  ratio  chi- 
square  test  of  the  fit. 

Converts  the  independent  variable  to  its  log  base  10  value  for  analysis. 
Computes  the  confidence  limits  for  the  independent  variable 


The  “class”  statement  on  the  next  line  names  the  response  variable  (mvl)  for  the  probit 
analysis.  The  “model”  statement  specifies  the  response  variable  and  the  independent 
variable  it  is  based  on.  The  “model”  statement  also  has  several  flags  associated  with  it  that 
specify  output  options  and  probit  parameters.  The  following  is  a  list  of  the  options 
specified  and  their  purpose. 


lackfit  Performs  a  Pearson  chi-square  test  of  the  fit  and  log-likelihood  ratio  chi- 

square  test  of  the  fit. 

d  Specifies  the  type  of  statistical  distribution  to  use. 

corrb  Prints  the  estimated  correlation  matrix  of  the  parameter  estimates, 

covb  Prints  the  estimated  covariance  matrix  of  the  parameter  estimates, 

inversecl  Computes  the  confidence  limits  for  the  independent  variable. 


The  last  statement  is  the  “ouQjut”  statement.  This  statement  associates  data  output  from 
the  probit  procedure  with  variable  names  that  can  be  used  later  in  the  SAS  run.  The  rest 
of  the  statements  in  the  script  specify  printing  options  for  the  graph  output. 
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APPENDIX  B: 

EZ-PROBIT  COMPUTER  PROGRAM 
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The  following  program  was  written  by  Captain  Lonnie  Manning  to  perform  the  probit 
calculations  as  ouflined  in  Finney’s  book. 

EZ-PROBIT 

TOPIC:  Overview  of  how  to  calculate  ED50  using  the  EZ-Probit  method. 

Starting  point  for  this  type  of  “statistic”  is  the  probability  function 

t  -  X 

P  (,  X  )  =  \  f  it )  d  t 


EZ-  Probit  assumes  a  normal  distribution  for  this  type  of  data.  Thus, 


fix)  = 


exp( 


2a 


Making  the  (probit)  transformations 


=  and  r  = 


we  get 


-n^Tc" 


The  normal  equivalent  deviate  (NED)  uses,  Y  =  ,  which  can  take  on  negative 

a 

values.  Y-5  typically  does  not. 


THE  CALCULATION 
THE  DATA 

Xi  Logic  of  the  i*"*  unique  energy  level. 
ni  The  number  of  trials  at  x, 
ri  The  number  of  lesions  for  Xi 

Pi  =  Estimated  Fix), 

m 
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WORKING  PROBITS 


The  equation,  7-5 


,  can  be  rewritten  as; 
a 


Yf  =a-X;  +b. 


The  basic  concept  behind  probit  is  to  find  the  slope  and  intercept  of  the  above  equation. 
Straight-forward  linear  regression  (LR)  or  weighted  LR  does  not  provide  the  most 
accurate  description  of  this  system.  This  is  accomplished  by  introducing  “working 
probits”  (WP). 

The  WP  is  derived  from  the  maximum  likelihood  equation  and  is  given  by: 


where 


(p,-P(Y,-5)) 


5 


1  .  (Y,-5)\ 

Zi  =-^exp( — 2^:^ — ). 


For  Pi  =  0,  we  have  the  famous 


yo-Yi- 


P(Yi-5) 

Zi 


For  Pi  =  1,  we  have  the  famous  yjoo 


3^100 


l-P(I^-5)  y.Q(Yi-5) 

’■  Z;  '■  Z,- 


PERFORMING  THE  CALCULATIONS 


Step  1. 

Choose  an  approximate  slope  and  intercept  for  the  line,  =aXi  +b. 
For  most  cases,  7^  =  5 ,  will  do  nicely  as  a  starting  point. 


Step  2. 

Calculate  a  new  slope,  a'  and  intercq)t  b'  using  weighted  LR  with  a  weighting 
factor  of: 


W;  = 


Z/ 


'  P(l^-5)(l-P(i;.-5)) 
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Step  3, 

Calculate  new  using  the  results  of  step  2. 

Yi'  =a'Xi+b\. 


Step  4. 


Do  steps  2  and  3  until  it  converges. 


Step  5. 


Once  the  slope  and  intercept  are  obtained.  Log  (EDso)  can  be  calculated. 
t=r-5  ,  1 

Since  P(Y  -  5)  =  f  -j=exp(-—)dt  =- ,  for  Y=5, 
jL  ^|2n  2  2 


5-b'  5-y+a''X 

Log(ED50)  =  XjQ  —  j  ,  or  x^q  —  j  . 

a  a 

Note:  In  reality,  one  calculates  the  equation,!"  =  y  +  -  x) ,  since  y,  a'  and  x  are  used 

later  to  calculate  FLs. 

THE  FOLLOWING  PROGRAM  PERFORMS  ALL  OF  THE  REQUIRED 
CALCULATIONS: 


Subroutine  1.  EZ-IOJl 

int  iofile  (long  double  *a,mt  *b, string  f_name  jnt  lin_log){ 
long  double  DUM_a  ,en; 

ifstream  inf(f_nanie.c_str(),  ios::in  I  ios::nocreate ); 

int  i4miaxj£jDUM_b,low4iit; 

if  (lin_log--l) 

{ 

for  (i-0;iiif.eof()!— l;i++) 

{ 

iiif»en; 
inf>>liit; 
if  (inf.eofO!“l) 

{  a[i]-logl01(en); 
b[i]-hit; 

> 


} 
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else 

{ 


} 


for  (i“0;mf.eof()!“l;i-H-) 

< 

inf>>en; 


inf>>hit; 

if(inf.eofO!-l) 

<  a[i]-en; 

b[i]-hit; 

> 


inf.closeO; 


iimax=i-2; 


for  (k'=0;k<-nmax;k++)  //Cheap  &  dirty  sorting  procedure 

{ 

low=k; 

for  (i“k+l;i<-nmax;i-H-) 

{ if(a[i]<=a[low])  low-  i;  }; 

DUM_a-a[low];DlJMJ)=b[low]; 

a[low]=aIk];b[low]-b[k]; 

a[k]-DXJM_a;b[k]=DUM_b; 

}; 

return  nmax; 

}; 


Subroutine  2.  N-INV.H 

routine  2.  //The  Inverse  of  die  standardized  Normal  distribution  use  for  fiducial  limits 


int  set_ninvGong  double  * ); 
int  set_ninvGong  double  *a  ){ 
a[0]-  -2.3263478740408411; 
a[l]-  -2.05374891063182305 
a[2]=  -1.88079360815125094 
a[3]=  -1.75068607125216998 
a[4]=  -1.64485362695147271 
a[5]-  -1.55477359459685354 
a[6]=  -1.47579102817917073 
a[7]-  -1.40507156030963255 
a[8]=  -1.34075503369021637 
a[9]=  -1.28155156554460047 


//norm  inverse  of  .01 
//norm  inverse  of  .03 
//norm  inverse  of  .05 


//norm  inverse  of  .10 


a[  1 0]=  - 1 .03643338949378955 ;  //norm  inverse  of .  15 


a[l  1]=  -0.841621233572914189; 

a[12]=  -0.674489750196081728; 

a[13]=  -0.524400512708040724; 

a[14]-  -0.385320466407567563; 

a[15]=  -0.253347103135799736; 

a[16]=  -0.125661346855073969; 

a[17]=  0.00000;  //norm  inverse  of  .5 

a[18]=  0.125661346855074109; 

a[19]-  0.25334710313579988; 
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a[20]-  0.385320466407567713; 

a[21]- 0.524400512708040884;  //by.05 

a[22]-  0.674489750196081877; 

a[23]-  0.84162123357291437; 

a[24]-  1.03643338949378979; 

a[25]-  1.28155156554460076;  //nonn  inverse  of  .90 

a[26]-  1.34075503369021653; 

a[27]-  1.40507156030963272;  //by  .01 

a[28]-  1.47579102817917092; 

a[29]=  1.55477359459685374; 

a[30]-  1.64485362695147295; 

a[31]=  1.75068607125217026; 

a[32]=  1.8807936081512513; 

a[33]=  2.05374891063182355; 

a[34]=  2.32634787404084201; 

return  0;}  ; 


Subroutine  3.  Probit«H 

#include  <iostreaniJi> 

#indude  <niathJh> 

#include  <stdioJi> 

#defme  PI  3.14159265358979323846264338327950 

//The  natural  log  of  the  gamma  function 

long  double  gamln(long  double  n) 

{long  double  sum; 

if(n— 0.5)  return  logl(sqrtl(PI)); 

if  (n— 1.0)  return  logl(l.O); 

if  (n=-2.)  return  logl(l.O); 

sum-0; 

if((n-(int)n)!-  0.5) 

{  for(int  i=0;i<— n-2.;i-H-) 

{ 

sum-sum+logl(n-l  .0-i*  1.0); 

} 

} 

else 

{ 

sum=logl(sqrtl(PI)); 
for(int  i-0;i<-(int)(n-1.0);i++) 
sum-sum+logl(n- 1 .0-i*  1 .0); 

} 

return  sum; 

} 

// _ _ 

//  The  normal  distribution  —  series  ^proximation 


long  double 
P_ser(long  double  x){ 
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long  double  af,an,sum,px,dj; 

af=0.;an-0.;sum=0.;px-0.;dj=0.; 

mtj^_x; 

m_x-l; 

if(x<0.0)  {x=-x;m_x=-l;> 

an-  1.; 
suni=an; 

for(j-l;j<-500;j++) 

{ 

af=sum; 

dj»*j*l.; 

an-an*x*x*(2.0*dj-l  .)/(dj*2.0+l  .0)/(2.*dj); 
sum=sum+an; 

if  ((sum-af)—!  .e-18)  break; };  /*  get  out*/ 

sum-x*sum/sqrtl(2.*PI) ; 

/*couLprecision(22);*/ 
px=l  ./2.+m_x*sum; 
return  px; 

} 

//Factorial  function  ~  only  good  for  integer  and  half-integer  values 

long  double 
fac(long  double  n){ 
if(n— 1.)  return  sqrt(-l); 
if  (n-=— .5)  return  sqrtlCPI); 
if  (n—O.)  return  1.; 
if  (n”l.)  return  1.; 
return  n*fac(n-l.); 

} 

//The  Gamma  function 
long  double 

gam(long  double  x){retum  fac((x-1.0));} 

// 

long  double 

gser(long  double  a, long  double  x) 

{ 

long  double  an,sum,eps; 
an-1.0; 
eps=l.e-20; 
sum-an; 

for(int  n=l;n<=5000;n-H-) 

{ 

an=x*an/(a+n*  1.0); 
sum=sum+an; 

if(fabsl(an)<fabsl(sum)*eps)  break; 

> 

return  sum*expl(-x)*powl(x,a)/a; 

} 
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// _ _ _ 

//calculates  p(chi,v)  using  the  continued  fraction  method 
long  double 

p_chi(long  double  CHI_2,int  df) 

{ 

long  double  aj,bj,cj,dj.fj.tiny,x,deltaJ,eps,ajes; 
intj; 

tiny=1.0e-45; 

eps-l.Oe-20; 

a-df/2.0; 

x-CHI_2/2.0; 

if  (x<a+1.0){res=  l.-gser(a^)/gam(a);goto  get_out;}; 
bj“X+1.0-a; 
if(bj==0.0)  bj=tiny; 
cj-l.O/tiny; 
dj=1.0/bj; 

forO‘“l  ;j<-5000;j-H-) 

{ 

bj+=2.0; 

dj=aj*dj+bj; 

if  (fabsl(dj)<tiny)  dj-tiny; 
cj-bj+aj/cj; 

if  (fabsl(cj)<tiny)  cj-tiny; 
dj-l.O/dj; 
deltaj=cj*dj ; 
fj=deltaj*fj; 

if(fabsl(delta_j-1.0)<eps)  break; 

} 

res=powl(x,a); 

res=res*expl(-x); 

//if(res<eps)  return  0.0; 

res=res/gam(a); 

rcs=res*fj; 

get_out: 

if  (res>1.0)  return  1.0; 
return  res; 

} 

// - - - 

//Probability  -Inf  to  x  for  the  Normal  distribution 

long  double 
PGong  double  x) 

< 

long  double  Qx; 

if  (fabsl(x)<=.6)  return  P_ser(x); 

Qx=  p_chi(x*x,l.)/2.; 
if  (x'^)  return  Qx; 
return  l.-Qx; 

} 
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//  The  inverse  of  the  Probability  for  the  Normal  distribution 
long  double 

norms_inv(  long  double  prob_x ) 

{ 

long  double  d,aO,al^aO; 

inti; 

aO-.l; 

for  (i=l;  i<=2000;i++) 

{ 

Za0-=expl(-a0*a0/2.0)/sqrtl(2.0*PI); 

al-aO-(P(aO)-prob_x)/ZaO; 

d=al-aO; 

if  (fabsl(d)<-l.e-18)  break; 
aO-al; 

}; 

//cout «  "d  ="«d«"  al>=  "«al«endl; 
return  al; 

} 

//  The  contnued  fraction  appraoximation  of  the  Incomplete  Bets  Function 
long  double 

betacf(long  double  a,long  double  b,long  double  x) 

{ 

int  m,Tn2; 

long  double  epsJFPMIN,aa,c,d,deld),qab,qam,qap; 
eps-l.e-20; 

FPMIN«l.e-45; 

qab=a+b; 

qap=a+1.0; 

qam=a-1.0; 

C”1.0; 

d=1.0-qab*x/qap; 
if(fabsl(d)<FPMIN)  d-FPMIN; 
d=1.0/d; 
h-d; 

for  (m-l;m<«=5000;m++) 

{ 

m2-2*m; 

aa-m*(b-m)*x/((qam+m2)*(a+m2)); 

d=1.0+aa*d; 

if(fabsl(d)^MlN)  d-FPMIN; 
c=1.0+aa/c; 

if(fabsl(c)<FPMIN)  c-FPMIN; 

d-l.O/d; 

h*=d*c; 

aa=-(a+m)*(qab+m)*x/((a+m2)*(qap+m2)); 

d=1.0+aa*d; 

if(fabsl(d)<FPMIN)  d-FPMIN; 
c=1.0+aa/c; 

if(fabsl(c)<FPMIN)  c-FPMIN; 

d=1.0/d; 

del-d*c; 
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h*-del; 

if  (fabsl(del-1.0)<eps)  break; 

} 

return  h; 

} 

//  The  incomplete  Beta  function 
long  double 

betai(long  double  a,  long  double  b,long  double  x) 

{ 

long  double  bt; 

if(x<0.0llx>1.0)  sqrtl(-l); 
if(x“0.0  II  x==1.0)  bt=0.0; 
else 

bt-expl(gamln(a+b)-ganiln(a)-gamln(b)+a*logl(x)+b*logl(1.0-x)); 

if(x<(a+1.0)/(a+b+2.0)) 

return  bt*betacf(a,b,x)/a; 
else 

return  1.0-bt*betacf(b,a,l'0"X)/b; 

> 

//  Prob  for  the  Students  t  distribution 
long  double 

student_t(long  double  t,long  double  v/*deg  of  freedom*/) 

{ 

return  1.0-betai(v/2.,lV2.,v/(v+t*t)); 

} 

//  Hie  inverse  of  Student's  t  distribution 
long  double 

student_t_inv(  long  double  prob_t,long  double  v ) 

{ 

long  double  d,a0,al^0,A00; 

inti; 

a0=0; 

A00=(0.5*logl(v)+gamln(0.5)+gamln(v/2.0)-gamln((1.0+v)/2.0)); 

//+((v+1.0)/2.0)*logl(1.0+a0*a0/v)); 

for  (i-1;  i<=500;i++) 

{ 

A0=A00+((v+l  .0)/2.0)*logl(l  .0+a0*a0/v); 

al-a0-(student_t(a0,v)-prob_t)*expl(A0)/2.0; 

if  (fabsl(al-a0)<=l.e-18)  break; 

aO-al; 

}; 

//cout «  "d  -"«d«"  al“  "«al«endl; 
return  al; 

} 
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Main  Program.  ProbitCPP 


#include  <ow^applicat.h> 

#iiiclude  <owl\framewin.h> 

#mclude  <owMialog  Ji> 

#include  <owNic.h> 

#include  <owr'opensave.h> 

#include  <owI'^Lh> 

#mclude  <owN:heckbox.h> 

#include  <owl\validate.h> 

#include  <c:Mx:4SmcludeV>wIVwindow.rii> 
#include  <stdioJi> 

#include  <stringdi> 
ffinclude  <ctype.h> 

#include  <math.h> 

#pragitia  hdrstop 

Ifeclude  "c:\bc4Smin.e'ez_ioJi" 

#include  "c:'45c4\imneNpobiLh" 

#include  "c:M3c4\mine\ii_mvJi" 


#define  CM_EMPINPUT  201 
#defineCM_FILEOPEN  0x100 
#define  CM.COLOR  0x101 

#define  CM_HELPABOUT  0x200 

long  double  PowL(long  double  xx,int  ln_lg){if  (ln_lg“=l)  return  powl(10.0pcx);retum  xx;}; 

long  double  ZQong  double  YY){ 

return  expl(-YY*YY/2.)/sqrtl(M_PI*2.); }; 

char  *fii; 

char  fln[24]; 

#define  MAXNAMELEN  35 
#define  MAXSSLEN  12 
#defme  MAXIDLEN  4 
#define  MAXDEPTLEN  7 
#defme  MAXSECLEN  3 

struct  TParamStruct  { 
char  DefaultDir[50]; 
char  Converg[MAXIDLEN]; 
char  OprecEdit[3]; 
char  Confidence[7]; 
char  PerCent[3]; 

BOOL  Log_Dose; 

BOOL  PennFile; 

BOOL  Exempt; 

}; 


// - 

class  TParamDlg :  public  TDialog  { 
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public: 

TParamDlgCrWindow*  parent,  const  char*  name,  TParamStruct&  tx); 

}; 


TParamDlg::TParamDlg(TWindow*  parent,  const  char*  name,  TParaniStruct&  tx) 
:  TDialog(parent,  name), 

TWindow(parent) 

{  TEdit  *edit; 

edit— new  TEdit(liiis,  101,  sizeof(tx.DefaultDir)); 

//edit->SetValidator(new  TFilterValidator("A-Za-z.'\N: ")); 
edit  -  new  TEdit(this,  103,  sizeof(tx.Converg)); 

edit->SetValida*<*r(new  TRangeValidator(12, 18)); 
edit  =  new  TEdit(this,  104,  sizeof(tx.OprecEdit)); 

edit->SetValidator(new  TRangeValidator(0, 21)); 
edit  =  new  TEdit(this,  105,  sizeof(tx.Confidence)); 

edit->SetValidator(newTPXPictureValidator("#.####")); 
edit  =  new  TEdit(this,  102,  sizeof(tx.PerCent)); 
edit->SetValidator(new  TRangeValidator(0, 100)); 
new  TCheckBox(1his,  106); 
new  TaieckBox(this,  107); 
new  TQieckBox(this,  109); 

//  TEdit*  edit; 

/*  edit  =  new  TEdit(this,  104,  sizeof(tx.DeptEdit)); 
edit->SetValidator(new  TPXPictureValidator("SalesJ)ev,Mfg"));  /* 
edit  —  new  TEdit(this,  105,  sizeof(tx.SecEdit)); 
edit->SetValidator(new  TPXPictureValidator("l  1,12,13,14,15")); 

*/ 

TransferBuffer=(void  far  *)&tx; 


} 

class  TCommDlgWnd :  public  TFrameWindow  { 
public: 

TCommDlgWndCTWindow*,  const  char* ); 
~TCommDlgWnd0; 

void  Paint(TDC&,  BOOL,  TRect&); 
void  CmFileOpenO; 
void  CmHelpAboutO; 
void  CmEmpInputO; 

TColor  Color; 

TFont*  Font; 

TOpenSaveDialog::TData  FilenameData; 
void  EvSize(UINT  sizeType,  TSize&  size); 
private; 

TParamStruct  ParamStruct; 

DECLARE_RESPONSE_TABLE(TCommDlgWnd); 

}; 
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DEFINE_RESPONSE_TABLEl(TCoimnDlgWnd,TFrameWindow) 

EV_WM_SIZE, 

EV_COMMAND(CM_FE^OPH^,  CmFileOpen), 

EV_C0MMAND(CM_EMPINPUT,  CmEmpInput), 
EV_COMMAND(CM_HELPABOUT,  CmHelpAbout), 
END_RESPONSE_TABLE; 


TCominDlgWnd::TConmiDlgWnd(TWmdow*  parent,  const  char*  title) 

:  TFrameWindow(parent,  title), 

TWindow(parent,  title), 

FilenaineData(OFN_FII^MUSTEXISTIOFN_HIDEREADONLYIOFN_PATHMUSTEXIST, 

"Data  Files  (*.dat)l*.datlAll  Files  (*  *)l*  *1", 

0, "*") 

{ 

memset(&ParamStruct,  0,  sizeof  ParamStruct); 


ParaniStructXog_Dose=TRUE; 

ParamStrucLPermFile-FALSE; 

strcpy(ParaniStrucLConverg,"18"); 

strcpy(ParamStructDefaultDir,"C:\^obit\v"); 
strcpy(ParamStruct.OprecEdit,'  '6"); 
strcpy(ParamStruct.Confidence,"0.1000"); 
strcpy(ParamStrucLPerCent,"95"); 
AssignMenu("CMDLGAPMENU"); 

//  Set  up  the  menu 

Color  =  TColor:  :Black;  //Use  black  as  the  default  color 

Font  -  0; 

} 


//  Empty  the  handle  to  the  font 


TCommDlgWnd:  :~TCommDlgWndO 

{ 

delete  Font; 

} 

// 

//  We  need  to  invalidate  the  entire  area,  not  just  the  clip  area  so  fliat 
//  paint  gets  called  correctly 

// 

void 

TCommDlgWnd:  :EvSize(UINT  sizeType,  TSize&  size) 

{ 

InvalidateO; 

TFrameWindow::EvSize(sizeType,size); 

} 


// 

//  Display  the  file  name  using  the  selected  font  in  the  selected  color. 

// 

void 

TCommDlgWnd::Paint(TDC&  paintDC,  BOOL,  TRect&) 

{ 

paintDC.SetTextColor(Color); 
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paintDC.SetBkColor(::GetSysColor(COLOR_WINDOW)); 
if  (Font) 

paintDC.SelectObject(*Font); 
paintDC.DrawText(fii,  strlen(fh), 

GetClientRectO,  DT_CENTER  I 

DT.WORDBREAK); 

} 

// 

// 

void 

TCoinniDlgWnd::CmFileOpenO 

< 

char  msg[50]; 

//MessageBox(itoa(ParaniStruct.Log_Dose,msg,10),"ASDASD",MB_OK); 

strcpyCmsg,"  "); 

fii-new  char  [512]; 

if  (TFileOpenDialog(this,  FilenameData).Execute()  -=  E)OK) 

{ 

strcpy(nisgJ’ilenameData.FileName); 

FlashWmdow(l); 

SetCursor(NULLJDC_WAIT); 

int  linjog,  PenaData; 

lin_log“ParamStrucLLog_Dose; 

long  double  confidence, tolerance, eps,outprec, percent; 

outprec=atoi(ParaniStruct.OprecE(iit); 

eps=powl01(-atoi(ParainStrucLConverg)); 
confidence-_atold(ParamStructConfidence); 
if  (confidence>1.0)  confidence=1.0; 
if  (confidence  <0)  confidence  -0; 

percent—atoiCParamStructPerCent)* 1 .0/100.00;  // outprec=  15; 

PennData-ParamStructPennFile; 

//change  to  default  to  6  for  clarence 
long  double  *  en;  en-  new  long  double  [200]; 

long  double  *  len;  len-  new  long  double  [2(X)]; 

long  double  *  y;  y-  new  long  double[200] ; 

//  long  double  *  Y;  Y=  new  long  double  [200]; 

long  double  *  p;  p=  new  long  double  [200]; 
long  double  *  ninv;  ninv=  new  long  double  [35]; 
int*  hit;  hit- new  int  [200]; 

int*n;  n=  new  int  [200]; 

int*H;  H- new  int  [200]; 

int*r;  r- new  int  [200]; 

long  double  *  w;  w-  new  long  double  [200]; 

long  double  YY,cr,sO,sl,s3,Sl,S2,S3,Sxx,Sxy,Syy,z,nw,x_,y_,chi_2xy,chi_2LR,gJT:.J),intercept; 

long  double  foo,oof,  slope,Yint,  EXP_P,yy,  log_ED50,  ED50,m,dum,tdiJ=LUfLL; 
int  j  JcJanax^imaXjbigslope^tot^ewt; 

newt=0; 

bigslope=0; 

string  filename,ofile,ofile2,ofile3,pnt,pth  Jidrl  4idr2; 
hdrl="Prob  Log[Dose]  Log[Lower  Limt]  Log[Upper  Limit]  V; 
hdr2-*Prob  Dose  Lower  Limit  Upper  Limit  V’; 
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set_ninv(ninv); 

pnt-""; 

filename-msg; 

strcpy(fii4ilename.c_strO); 

for  (k-0;k<-filename.ifed('\\");k++)  pth=  pth+filename[k]; 
strcpy(ParamStrucLDefaulfDir,pth.c_strO); 
for  (ic-0;k<»filename.rfind(".");k-H-)  pnt-  pnt+filename[k]; 

if(PeniiData=FALSE)  pnt="tmp."; 
ofile-pnt+"out"; 


if(lin_log=l) 

<  ofi]e2-pnt+'log"; 

ofile3-pnt+"grf'; 


} 

else 

{hdrl=hdr2; 

ofile2-pnt+"grf"; 

} 

ofstream  otf(  ofile.c_strO,  ios::out ); 
ofstream  otf2(  ofile2.c_str0,  ios::out ); 
ofstream  otf3(  ofile3.c_str0,  ios::out ); 
otf.precision(outprec); 

//  otf«ofile  «Nn'; 

otf.precision(outprec); 

nmax-iofile(en4iit^ename,lin_log); 


/*the  following  procedure  acts  like  a  "pivot  table"*/ 
r[0]-hit[0];  H[0]-0;  k-1;  foo-en[0];  ntot-hit[0]; 
for  (i-l;j<-nmax;j-H-) 

{  ntot-ntot+hitlj]; 

if(en[j]“foo) 

r[k-l]=hit[j]+r[k-l]; 
else  /*  new  energy*/ 

{ 

H[k]-j; 

r[k]=liit[j]; 

k++; 

foo-en[j]; 

} 

} 

kmax-k-1; 

H[kmax+l]=nmax+l;  /*need  to  make  sure  H[kmax+1]  exists*/ 
for  (k-0;k<-kmax;k-H-) 

{ 

n[k]-H[k+l]-H[k]; 

j-H[k]; 

len[k]-en|j]; 

} 

//Initial  values  for  slope  and  Intercept  of  the  Probit 
//for  normit  Yint  —0 
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slope=0.0; 

Yint=5.0; 

log_ED50=10.0;  //log_ED50-ED50  for  lin_log=TAKE_LIN 

dum-1.0; 

t=nonns_inv((long  double)0.975); 
do{ 

oof-log_ED50; 

s0“0;sl“0;s3»^;Sl— 0;S2“0;S3“0,chi_2LR“0; 
for(k=0;k<=kmax;k-H-) 

{ 

p[k]-(1.0*r[k])/(1.0*n[k]); 

yy-lenDc]*slope  +Yint;  //yy-expected  probit 

EXP_P=P(yy-5.); 

YY=yy-5.0; 

foo“p[k]-EXP_P; 

if(foo==0.0) 

{  yM=yy;  //lim  as  YY->INF  foo/z  =0, 2f=0 
w[k]-0.0; 

} 

else 

< 

YY-YY*YY/2.0; 

Z“expl(YY) ; 
z-z*sqrtl(M_PI*2.0); 
z»1.0/z; 
yM-yy+foo/z; 

w[k]-z*z/(EXP_P-EXP_P*EXP_P); 

} 

nw=n[k]*w[k]; 

sl-sl+nw*len[k]; 

51- Sl+nw*lenDc]*len[k]; 
s3=s3+nw*y[k]; 

52- S2+nw*y[k]*len[k]; 

53- S3+nw*y  [k]*y[k] ; 
s0=s0t-nw; 

//dii_2LR=chi_2LR+(r[k]*1.0-  n[k]*EXP_P)*(r[k]*l.a- 
n[k]*EXP_P)/n[k]/(EXP_P-EXP_P*EXP_P); 

}; 

Sxx-Sl-sl*sl/sO; 

Sxy-S2-sl*s3/s0; 

Syy=S3-s3*s3/s0; 

slope-Sxy/Sxx; 

x_-sl/s0;  //the  new  equation  is  Y-y_+slope*(x-xJ; 

y_=s3/s0; 

Yint-y_-slope*x_;  //i.e.  X“0 
log_ED50=x_+(5.0-yJ/slope; 
if(slope>200.)  {bigslope“l;break;} 
if(dunr>55)  break;  //Maximum  number  of  iterations 

dum-dum+1.0; 

//tolerance  =  fabsl(oof-oof); 
tolerance  =  fabsl(Gog_ED50-oof)); 

}  while  (tolerance>-eps); 
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chi_2xy=Syy-Sxy*Sxy/Sxx;  //CHI  Squared 
in=log_ED50;  //logEDSO  or  ED50 

if(chi_2xy<l  .Oe-18)  //Chi'^2  -  zero 

{chi_2xy»=0.0;oof- 1 .;  > 
else 

oof=p_chi(chi_2xydanax-l);  //Probability  of  Chi'^2,  N-2  Degress  ofFreedom 
h-1.0; 

if(fabsl(percent— 0.95) ) 

t=1.96;  //t-student_t_inv(percent,INF); 

else 

t=nornis_inv(l.-(l.-perceiit)/2.);  //Calculates  t  for  new  confidence  level 

if(oof<confidence)  //Check  for  good  data 

{newt=l; 

h=chi_2xy/(kmax*  1 .0- 1 .0);  //heterogeneity  factor 

t-student_t_mv(percent,(long  double)(lanax-l)); }; 
intercept=y_  -slope*x_-5.0; 

ED50=PowL(log_ED50,lin_log); 

//cr=Sxy/sqrtl(Sxx*Syy) ; 

//cr-cr*l.ci0; 

oops: 

if(bigslope!-l)  //Fix  for  infinite  slope 

< 

g=(t’'‘t)*h/Sxx/slope/slope; 

FL=m+g*(m-x_)/(1.0-g); 

D-  (l.-g)/sO+(ni-x_)*(m-x_)/Sxx; 
if(D<0.0){bigslope-l;goto  oops;} 

D-sqrtl(h)*t/slope/(1.0-g)*sqrtl(D); 

FLL”FL-D; 

FLU-FL+D;> 

else 

{h=0; 

g=0; 

FLL=m; 

FLU-m;} 

//I  I  I  I  I ++++++++++M  ++-I-++++++++-H-H  I  I  I  I  I  I  n  I  I  I  I- 1  >+++ 1  1  + I  I  t+  l  I- »++++++++++++++++ 
++ 

otf«"  "«filename.c_strO«"\n\n" 

«'Nt\t  ONES  =  "«ntot«'Nfi  ZEROES  -  "«(mnax+l-ntot)«V 
«'M'^  h-"«h«V 

«'Nt\t  g  -"«g«'\n' 

«'\fst  t  -"«t«V 

/*«'Nt\t  Snw  =  "«s0«^o' 

«"\t\t  Snwx  -  "«sl«^n' 

«'\t\t  Snwy  -  "«s3«Nn' 

«'NtSt  Snwxx-"«S1«V 
«'\t\t  Snwxy  ■=  "«S2«V 

«'Nt\t  Snwyy  -  "«S3«^o'*/ 

«'Nt\t  SYY  -  "«Syy«V 

«'Nt\t  SXY  =  "«Sxy«V 

«’'^\t  SXX-"«Sxx«V; 

if(lin_log“l)  otf«'Nt\t  LGED50  -  "«log_ED50«Nn' ; 
otf«'M\t  ED50-"<<ED50«V 

«'Nt\t  FLU  -  "«PowL(FLU,linJog)«V 
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«'Nt\t  FLL  -  "«PowL(FLL,lin_log)«"^nV 

«"  Pearson's  Cffl2  -  "«chi_2xy«"St  DF  -  "«(kmax-l)«V 
«"  Prob>Chi-Sq  -  "«  oof«"yriSn" 

«"M\t  XBAR  =  "«x_«'Nn'  //only  kmax-l=N-2  because 
YBAR  -  "«y_«'Sn'  //  k  goes  from  0  to  kmax 

«'NtNt  Intercept  -  "«mtercept«V 
«"\bt  Slope  -  "«slope«'\n' 

«'lteration  -  "«(mt)  dum«"Sn\n\nNn"«endl; 

//+■!  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  M  I  n  -i  1 1  I  ++++I  I  I  H  +++  I  I  I  I  I  I  1  +  I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  1+++ 
+++ 

otf2«'  Nn  Input  File  "«filename.c_str()«"\n\n" ; 

otf2«hdrl  ; 

//  +0.01  -0.881940268668475  -1.78261681790058  -0.572033298124713 

otf2.precision(outprec); 

otf2.setf(ios::showposlios::left); 

otf2«'\n'; 

for(k- 1  ;k<- 1 0;k++) 

{  m-(ninv[k-l]-mtercept)/slope; 

FL-m+g*(m-xJ/(l.-g); 

D-t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 
otf2.width(5); 
otf2«(k*0.01)«V; 
otf2.widfli(18); 
otf2«m«’\f; 
if(bigslope!=l){ 
otf2.widlh(18); 
otf2«(FL-D)«Nt'; 
otf2.width(18); 
otf2«(FL+D)«'^n'; } 
else 
{ 

otf2.width(18); 
otf2«"  -  "«'St'; 
otf2.widtii(18); 
otf2«"-"«Nn';} 

}; 

for(k- 1  ;k<-l  6;k++) 

{  m-(ninv[k+9]-intercept)/slope; 

FL-m+g*(m-xJ/(l.-g); 

D-t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 

otf2.width(5); 

otf2«(.l+k*0.05)«V; 

otf2.width(18); 

otf2«m«'^t'; 

if(bigslope!=l){ 

ot£2.width(18); 
otf2«(FL-D)«'^; 
otf2.width(18); 
otf2«(FL+D)«Nn'; } 
else 
{ 

otf2.width(18); 

otf2«"-"«V; 
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otf2.widfli(18); 
otf2«"  - 


}; 


for(k-l  ;k<-9;k-H-) 

{  m=(mnvlk+25]-intercept)/slope; 

FL-m+g*(m-x J/(l  .-g); 

D=t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 
otf2.width(5); 
otf2«(.9+k*0.01)«V; 
otf2.width(18); 
otf2«m«V; 
if(bigslope!-l){ 
olf2.widlli(18); 
otf2«(FL-D)«Nt'; 
otf2.width(18); 
otf2«(FL+D)«'Sn'; } 
else 
< 

otf2.widfli(18); 
otf2«"  -  "«V; 
otf2.width(18); 

}; 

// - 

if(lin_log--l)  //do  this  part  because  we  also  want  die  file  in  regular  non-log  units 

{ 

//■l  ■l■++++^  I  I  I  I  I  <■  h++++++Ouqnit  tof  ile  *fh.grf  M +  i  1 1  ill  t  -n-++++++-i' <■++++++ 1■^■++++ 

otf3.precision(outprec); 

otf3«"^n  Input  File  "«filename.c_str()«'Nn'Sn" ; 

otf3«hdr2 ; 

//  +0.01  -0.881940268668475  -1.78261681790058  -0.572033298124713 

otf3.setf(ios:;showposlios:;left); 
for(k-l  ;k<-10dt++) 

{  m=(ninv[k-l]-mtercept)/slope; 

FL-m+g*(m-xJ/(l.-g); 

D-t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 

otf3.width(5); 

otf3«(k*0.01)«V; 

otf3.width(18); 

otf3<<PowL(m,lin_log)«'\t'; 

if(bigslope!“l){ 

otf3.width(18); 

otf3«PowL(^-D),lin_log)«V; 

otf3.width(18); 

otf3«PowL((FL+D),lin_log)«  V; } 
else 
{ 

otf3.width(18); 
otf3«’'  -  "«M'; 
otf3.width(18); 
otf3«"  -  "«^\n';} 
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}; 


for(k=  1  ;k<=  1 6;k++) 

{  m-(ninv[k+9]-mtercept)/slope; 

FL=m+g*(m-xJ/(l.-g); 

D-t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 

otf3.width(5); 

otf3«(.l+k*0.05)«V; 

otf3.width(18); 

otf3«PowL(m,lin_log)«V; 

if(bigslope!“l){ 

otf3.width(18); 

otf3«PowL((FL-D),liii_log)«V; 

otf3.widfli(18); 

otf3<'<PowL((FL+D),lin_log)« V; } 
else 
{ 

otf3.width(18); 

otf3«"-"«M'; 

otf3.width(18); 

otf3«"-"«V;} 

}; 


for(k=l;k<=9;k-H-) 

{  m-(ninv[k+25]-mtercept)/slope; 

FL-m+g*(m-x J/(l  .-g); 

D-t/slope/(l.-g)*sqrtl((l.-g)/sO+(m-xJ*(m-xJ/Sxx); 

otf3.width(5); 

otf3«(.9+k*0.01)«V; 

otf3.width(18); 

otf3«PowL(m,lin_log)«Nt'; 

if(bigslope!-l){ 

otf3.widtii(18); 

otf3«fowL((FL-D),lin_log)«V; 

otB.width(18); 

otf3«PowL((k.+D),lin_log)«V; } 
else 
{ 

otf3.width(18); 

otf3«"-"«V; 

otf3.width(18); 

otf3«"-"«V;} 

}; 

}; 

// - 

//.|.  I  .H-+4  I  I  I  I  I  I  M  I  I  >++Out  put  to  Main  Window  l■■^++^-++^  i  i  i  i  ( ++++++++++++-I  I  l  i  I +++++++ 

strcat(fii,"Nn  Intercept  -  ")  ; 

strcat(fn,gcvt(intercept,21,msg)); 

strcat(fn,'Nn"); 

strcat(fn,'\a  Slope  -  ")  ; 
strcat(fn,gcvt(slope,21  ,msg)); 
strcat(fQ,'Nn"); 
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{  strcat(fii,"^  LogED50  =  ") ; 

strcat(fh,gcvt(log_ED50,21  ^g)); 
strcat(fii,"vii"); 

}; 

strcat(fn,'Nn  ED50  =  ") ; 

strcat(£ii,gcvt(ED50,2l4iisg)); 
strcat(fh,'Nn"); 

strcat(fii,'Nn  Log  Fiducial  Limits  ") ; 
strcat(fD,'V'); 

strcat(fii,gcvt(FLL,2l4iisg)); 
strcat(fh,"  "); 
strcat(fh,gcvt(FLU,21  ,msg)); 
strcat(fh,'Nn"); 

otf3.close(); 

otf2.close0; 

otfxloseO; 

FlashWindow(0); 

InvalidateO; 

SetChirsor(NULLJDC_ARROW); 

ofile="notepad  "+ofile; 
ofile2-"notepad  "+ofile2; 

ofile3="notepad  "+ofile3; 


//  ofile2-"c:'^ExceI\'^xcel  c:\«xcer^robitxlt";  //this  works 
WiiLExec((coiistchar*)ofile.c_str(),SW_RESTORE); 

WiiiExec((const  char*)ofile2.c_strO,SW_RESTORE); 

if(lin_log--l)  Wii]Exec((const  char*)ofile3.c_str(),SW_RESTORE); 


} 

// 

// 

void 

TCoimnDlgWnd::CmHelpAboutO 

{ 

MessageBox("  Finney's  Probit  NnWritten  by  Dr.  Lonnie  W.  ManningV 

"  Copyright  (c)  1995  Nn  The  NutHouse", 
"About  EZ-PROBIT", 

MB_OK); 

} 

void 

TCommDlgWnd-CmEmpInputO 

{ 

char  empInfo[sizeof(TParamStruct)+2+l]; 
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if  (TParamDlg(diis,TARAMINFO",  ParamStruct)^xecuteO  ==  IDOK)  { 

strcpy(empInfoJ*araniStruct.DefaultDir); 
strcat(empInfo,  'V'); 
strcat(empInfo,  ParamStnicLConverg); 
strcat(empInfo,  'Nn"); 

/*  strcat(empInfo,  ParamStructDeptEdit); 
strcat(einpInfo,  "Nn"); 
strcatfemplnfo,  ParamStructSecEdit); 
strcat(empInfo,  'V);  *! 

H  strcpy(enipInfo,  ParanaStructXog_Dose  ?  'EuUTime  " :  "PartTime  "); 
strcat(empInfo,  ParaniStructXog_Dose  ?  "Log  Dose\tt  "  :  '1.111  Dose^n  "); 

strcat(empIiifo,  ParamStrucLPennFile  ?  "Save  OutputV :  'DonVt  Save  OutputNn  "); 
strcat(empInfo,  ParamStruct^xempt  ?  'Exempts  " :  "NonExemptNn  "); 
MessageBoxfempInfo,  'Information  Stored",  MB_OK); 

} 

} 

// - - - 

class  TComniDlgApp :  public  TApplication  { 
public: 

TCommDlgAppO :  TApplicationO  {> 
void  InitMainWindowO  { 

MainWindow  -  new  TCommDlgWnd(0,  'EZ-PROBIT"); 
MainWindow->SetIcon(this,  'EZ_PROBIT"); 

EnableCtlSdO; 

} 

>: 


int 

OwlMain(int  /♦argc*/,  char*  /*argv*/  []) 

{ 

return  TComniDlgApp().Run(); 

} 
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